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Abstract 


\ 


This report presents the results of a study of the application of and /i synthesis techniques 
to the design of feedback control laws for the longitudinal dynamics of the High Angle of 
Attack Research Vehicle (HARV). The objective of this study is to develop methods for 
the design of feedback control laws which cause the closed loop longitudinal dynamics of 
the HARV to meet handling quality specifications over the entire flight envelope. Control 
law designs are based on models of the HARV linearized at various flight conditions. The 
control laws are evaluated by both linear and nonlinear simulations of typical maneuvers. 
The fixed gain control laws resulting from both the H 0 0 and fi synthesis techniques result in 
excellent performance even when the aircraft performs maneuvers in which the system states 
vary significantly from their equilibrium design values. Both the H ^ and (i synthesis control 
laws result in performance which compares favorably with an existing baseline longitudinal 
control law. 


1 Introduction 

The overall objective of this research program is to develop control law design techniques for 
aircraft which operate in regimes in which nonlinearities in the equations of motion have an 
important effect on the dynamic response characteristics. Our long term goal is to develop a 
design methodology which combines nonlinear dynamic inversion and fi synthesis techniques. 
This should result in nonlinear control laws with performance and stability characteristics 
which are robust to modeling errors and parameter variations. For the initial stage of this 
research, we studied the application of H «, and fi synthesis techniques to the design of 
longitudinal control laws for the HARV. This report summarizes the results of this study. 

The open loop longitudinal dynamics of the HARV do not meet standard aircraft handling 
quality requirements at low speeds and relatively high angles of attack. Our design objec- 
tive is to develop control laws which result in closed loop dynamic response characteristics 
which meet these requirements. These control laws must not result in saturation of actuator 
deflection and deflection rates. The /r synthesis control laws are explicitly designed to be 
insensitive to (1) parameter uncertainties in the linear design models, (2) uncertainty in 
the effectiveness of the actuators, (3) uncertainty in the sensor measurements and (4) input 
uncertainties including higher order dynamics. 

Our basic philosophy is to design control laws which minimize the weighted norm of the 
difference between the actual aircraft pitch rate, angle of attack, and/or normal acceleration 
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responses to pilot stick inputs and the desired responses to these inputs as given by models 
based on standard handling quality specifications. Furthermore we wish to perform this 
minimization without saturating our control actuators. To achieve this goal we minimize 
the weighted H norm of the actuator deflection and deflection rates. 

We learned two important lessons which have general applicability to H <*> based techniques 
for aircraft control law design. First, we found that it was important to include phugoid 
dynamics in our model of desired aircraft handling qualities. H based techniques, including 
fj, synthesis, minimize the norm of the weighted system output over the entire frequency 
range. In the case of aircraft handling qualities enhancement, part of the output vector is the 
difference between some of the actual aircraft states and the desired values of these states as 
predicted by models which reflect the handling quality specifications. If there is a significant 
difference between the desired and actual aircraft response in a particular frequency range, 
Hoo based techniques will produce control laws which act to minimize this difference. If 
this difference occurs in a frequency range in which we are not interested in optimizing 
performance, the resulting control laws may not give very good results in frequency ranges 
in which performance is important. Most aircraft handling quality specifications are given in 
terms of short period response characteristics. If we formulate these characteristics in terms 
of transfer functions, the low frequency characteristics will in general be different from those 
of the actual aircraft since the phugoid mode has been neglected. In the case of angle of 
attack response, the phugoid mode cancels out of the transfer function and the low frequency 
response from the handling quality model and the aircraft axe essentially identical. This is 
not the case for the pitch and normal acceleration responses and the low frequency gains 
for the aircraft are considerably less than those of the short period handling quality models. 
In order to achieve good tracking of pilot inputs, it is essential to include some model of 
the phugoid dynamics in the handling quality models. It might be possible to achieve good 
tracking by suitable choice of output error weighting functions, i.e. reduced error weighting 
at low frequency, but the addition of a phugoid mode to the handling quality model is a 
natural way of addressing the problem which allows us to achieve both short period and 
phugoid handling quality specifications. This can be important at flight conditions at which 
the phugoid is lightly damped or unstable. 

Second, we found that is important to add a pilot stick shaping or pre-filter following the 
pilot input. Otherwise the control law acts to provide performance at high frequencies at 
which precise tracking is not needed at the expense of performance in the low frequency 
range in which tracking is important. 

Although the design techniques described in this report can be applied to a variety of handling 
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quality objectives, we concentrate on having the pitch rate follow pilot stick inputs at low 
angles of attack and the angle of attack follow pilot inputs at high angles of attack. 
and ji synthesis control laws are designed for linearized models of the HARV at several 
flight conditions. These control laws are evaluated using linear and nonlinear simulations 
of the HARV. The nonlinear simulations include lateral dynamics and the baseline lateral 
control laws included in the simulation are used. Both the and // synthesis control laws 
follow pilot inputs well and meet handling quality specifications. In the linear simulations 
the \i synthesis control laws exhibit less sensitivity to parameter variations than do the 
control laws, however, both sets of control laws exhibit almost identical responses in nonlinear 
simulations in which transient values of the state vary substantially from their equilibrium 
values. Both sets of control laws yield responses which are similar to those obtained with 
the baseline longitudinal control laws which were included with the simulation. 

The remainder of this report is divided into six sections. The next section gives a brief 
overview of the H and fi design methodologies. The third section describes the handling 
quality models used. These are based on recent handling quality specifications [1], The fourth 
section details the H oo design and includes a method for optimal allocation of control effort 
between the stabilizers and the thrust vector control (TVC). The fifth section describes the /i 
synthesis design. The sixth section introduces a comparison between the performance of one 
of our control laws and the baseline controller. The seventh section consists of conclusions 
and future research directions. 


2 Controller Synthesis Techniques 

2.1 Hoo Controller Synthesis 

In recent years there has been great interest in controller design based on Hoo control theory 
[2-7]. A substantial contribution to the solution of the Hoo problem has resulted from the 
state space representation of the problem [5]. The Hoo space is a Hilbert space [8] consisting 
of the bounded and stable functions. The Hoo norm of a real rational matrix T(s) is defined 
as: 


Ill’ll. = Sup ||T(,)|| = Sup ||T(ju)|| = Sup *[Ttf«0] (1) 

Rc(*)> 0 W W V 7 

where a[T(ju)] denotes the largest singular value of T(ju}) over frequency. Restricting the 
plants to Linear Time Invariant (LTI), the basic system, called P, considered in the H^ 


3 


synthesis is described by : 


( 2 ) 

(3) 

(4) 


x(t) = A x(t) + B\w{t) + B2u(t ) 
z{t) = C\x(t) + + D\2u(t ) 

y{t) = C 2 x(t) + D 2 iw(t ) + D 2 2ii(t) 

where x(t ) is an n dimensional state vector, w(t ) is an m\ dimensional disturbance vector, 
u(t ) is an m2 dimensional control vector, z(t ) is a p\ dimensional output vector whose 
magnitude is to be minimized and y{t ) is a P2 dimensional measurement vector. 

The objective of Hoo control synthesis is to find a stabilizing controller K, whose input is y(t) 
and whose output is u(t), that will minimize the Hoo norm of the transfer function between 
w(t) and z(t). When such a controller is found, the following condition is satisfied : 

Moo <7 (5) 

with 7 being the least upper bound of the Hoo norm. Solution for the Hoo controller involves 
the solution of two Algebraic Ricatti Equations (A.R.E) [4,5]. The first A.R.E to be solved is 
related to the full information case represented by the Hamiltonian matrix loo • The second 
A.R.E is related to the output estimation case. This A.R.E is represented by the Hamiltonian 
matrix Joo- The matrices /<» and Joo are given below. 

[ A T 2 BiB' x -B 2 B' 2 
00 " [ -C[Ci -A' 

[ a' r 2 c[Ci - c' 2 c 2 

00 [ -BiB[ -A 

The software used for controller synthesis [2] implements the bisection method on 7 to find 
the admissible controller which simultaneously satisfies equation (5) and the following three 
conditions: 

• loo £ dom(Ric) and Xoo = Ric(Ioo) > 0. 

• Joo € dom(Ric) and Yoo = Ric(Joo) > 0. 

• p(XooYoo) < 7 2 

Where p stands for the spectral radius of a matrix. A matrix belongs to dom(Ric) if 1) None 
of its eigenvalues are on the imaginary axis and 2) it has an n dimensional invariant 
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subspace X-H which corresponds to the n L.H.P eigenvalues which can be partitioned to 
X-H = (Xi X 2 )' with Xi being invertible. 

In such a case, the controller is defined by the following : 


K(s) 


' Aoo 

ZooL^ 

Foo 

0 


( 8 ) 


where: 

Am = A + 7 2 BiB 1 X 00 + B 2 F^ + ZoqZ-oo^ 
Foo - -B' 2 X oo 

Loo = -VooC, 

Zoo = (/ - T-^oo^oo )- 1 


In addition to the constraints given above, the system must meet the following 4 conditions: 


1 . 

2 . 

3. 

4. 


(A,B 2 ) is stabilizable and (C 2 ,A) is detectable. 
D \ 2 and D 2 \ have full rank. 

A — jul B 2 


C x 


D 
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A — juil B\ 
C 2 D 2 \ 


has full rank for all u. 


has full rank for all u. 


The first assumption is necessary for the existence of a stabilizing controller. The second 
means that the penalty on z includes a nonsingular normalized penalty on the control u, w 
includes both plant disturbance and sensor noise, and the sensor noise weighting is normalized 
and nonsingular. 

The third and fourth asumptions together with the first one guarantee that the two Hamil- 
tonian matrices involved in the solution belong to dom(Ric). These two assumptions are 
necessary for the solution described above. 

The assumption of unstructured disturbances results in the Hqo synthesis producing a con- 
troller which is usually conservative in the sense that the performance and robustness of the 
closed loop system are underestimated. The more realistic case of structured uncertainties 
is taken into consideration with // analysis. The scaling matrices from /x analysis make it 
possible for the designer synthesize a less conservative controller than is obtainable from the 
standard Hqo formulation. The use of /x for analysis and then for synthesis is described in 
detail in the following two sections. 
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Figure 1: The Linear Fractional Transformation (LFT) M = T zw = Fi(P, K) 

2.2 \x Analysis 

Let us first define p for a given matrix M € C n x n . In the definition of p(M) there is an 
underlying structure A £ A where A is a prescribed set of block diagonal matrices. We will 
partition A into the following form: 


A = { diag(8il r i , ,6 s I r3 ,Ai, , A^) : <$. G C, Aj E C 


m , X m, 


} 


( 10 ) 


where r stands for the repeated scalar blocks, A for a full matrix block, S is the number of 
repeated scalar blocks and F is the number of the full blocks. S and F satisfy the relation 
Ef=i r i + i m j — n i with n being the dimension of M as defined above. 


For the above M and the prescribed A, ^a(M) is defined [9, 10] by: 
fi A (M) = [min ct( A) : A € A, Det(I - MA) = 0 ] _1 


( 11 ) 


If no A makes (/ — MA) singular then fi A (M) = 0. It is shown in [9] that the following 
relations hold : 

max p(QM ) < max p(AM) = p A (M) < * n / a(DMD~ l ) fto) 

QeQ AgA d € D ' / 


where 

Q = 

{Qe A : 

Q'Q = In] 


D = 

{diag[ A u . 

• ••) Aj, ..., dp'Ifnp ] : 

(13) 


Di e C T ' x 

r * , Di = D* > 0 , dj 6 R , dj > 0} 
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After applying a feedback controller K to the system P, defined in equations (2)-(4), as 
shown in Figure 1, to get M = F\ ( P, K ) , we define a Linear Fractional Transformation 
(LFT) [3, 4] of M with the blocks Ai and A 2 . The matrix M can be partitioned as: 


M = 


Mn M\2 

M 2 i M 22 


(14) 


Where Mu is related to the parametric and model uncertainties and M 22 is related to the 
system performance. The dimensions of Ai and A 2 are compatible with the dimensions of 
Mu and M 22 respectively. A third structure is defined as: 


A = { 


A t 0 

0 A 2 


Ai e Ai,a 2 <e a 2 } 


(15) 


By the small gain theorem the LFT Fi(M, A 2 ) is well posed if and only if (/ — M 22 A 2 ) is 
invertible, where 


Fj(M, A 2 ) = Mu -b Mi 2 A 2 (7 — M 22 A 2 ) 1 M 2 i (lb) 

Recalling the definition in equation (11) for p, it can be seen that the LFT -F;(M, A 2 ) is well 
posed for all A 2 € A 2 if and only if /i 2 (M 22 ) < 1. // 2 denotes the fi test with respect to the 
block structure A 2 . 

The fi main loop theorem [3, 4, 9] says that the following are equivalent: 


Ha (M) < 1 /i 2 (M 22 ) < 1 and {fii(Fi(M, A 2 )) < 1} (17) 

Za 2 G £ 

The fi test defined at the left is the one that captures the robust performance of the sy- 
stem. The software we use performs an equivalent frequency domain fi test [2]. This test 
enables us to analyze the robustness of a given system in the presence of a prescribed set of 
perturbations. The fi analysis results make it possible, by scaling the synthesized system, 
to improve the results from the Hoc synthesis by making 7 in equation (5) smaller. This 
technique of scaling the system and then synthesizing an Hoc controller for the scaled system 
is the synthesis technique which will be described in the following section. 
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2.3 // Synthesis. 

The approach taken to fi synthesis in this paper is called D - K iteration. In D - K iteration 
the peak value of of the closed loop transfer function from w to z is minimized over 

all admissible controllers, i.e: 


7TIITI 

stabilizing {||F;(P, A')|| Ma } 
K(s) 


(18) 


where 


||G|U {a*a[G(jw)]} 


(19) 


We replace (Af) by its upper bound a(DM D~^) and assume that D € D is the appropriate 
scaling for the perturbed system. Recalling the upper bound in (12) and reformulating the 
fj, synthesis objective in (18), we get 


min { max { min a[D u Fi(P,K)^j w )D ul 1 ]}} 

Stabiliztnp uj D w gD 


( 20 ) 


where the inner minimization is the approximation to [Fi(P,K)p u )\. D w is the D sca- 
ling matrix which corresponds to Fi(P, K)^^) at the point u. Since D w is chosen from D 
independently at every frequency, equation (20) can be rewritten as: 


min { min { max a[DuFi(P, A/]}} 

Stabilizing u> 

K(s) 


( 21 ) 


or 


min { min \\D w Fi(P, ||oo} 

Stabilizing 

K(s) 


( 22 ) 


Consider a real-rational, stable, minimum-phase D(s) to approximate the magnitude of D w 
over frequency. So, (22) becomes: 


min { min \\D(s)Fi(P, K)fj^)D(s) 1 <»} 

Stabilizing £>{.) 

K(.) £>(.)eD 


( 23 ) 


The D-K iteration procedure is defined in (23) and the steps involved are: 
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1. Define a stable non minimum phase function D(s ) to fit the magnitude of D w 
that was computed through the n analysis. 

2. Use fifoo synthesis to find a controller K(s) which minimizes the scaled system 
norm ||Z)(s)F/(P, A')D(s) _1 || 00 . 

3. Run the jx test. Check n\Fi(P,K)] and build a new D w . 

4. Go through steps 1 through 3 until no further reduction in /x is achieved. 


3 Handling Qualities Requirements 


The Hqq and /x techniques require the characterization of the desired handling qualities of 
the aircraft to be defined in the frequency domain. Our objective is for the aircraft to satisfy 
level 1 handling qualities as defined in Ref. [1]. 

In defining the handling qualities we will use the short period linearized approximation of 
the aircraft dynamics given by: 

q = M q q + M a Q + MsS ^4) 

a = q + Z a a 


Where q is the pitch rate, a is the angle of attack, 8 is the stabilizer deflection and M q , M a , Ms 
and Z a are the aerodynamic derivatives. The first three are the moment derivatives with 
respect to pitch rate, angle of attack and stabilizer deflection respectively. The fourth aero- 
dynamic coefficient is the Z force derivative with respect to the angle of attack. The Z force 
due to the stabilizer deflection is neglected. 

The normal acceleration at the center of gravity, in the absence of lateral velocity, roll and 
yaw rates, is given by: 


Uq — w 

Tlzcg = 

9 


(25) 


where U and w denote the longitudinal and the vertical velocity components of the aircraft 
in body axes. Using (24) and (25) we get the steady state ratio of n Zcg to a as: 


( 


n zcg 

a 


)s$ 


UZg 

9 


with n z given in g units. 


( 26 ) 
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Assuming to be constant, the transfer functions from the stick to q, n Zcg and a are: 


stick 


n, 


~cg 


= K, 


stick 


1 _ i 

q i + + 4 


v *p 


w 


sp 


= K, 


nz 


1+^+4 


(27) 


(28) 


a 


stick 


= I<a 


l + ^2l + jL 


(29) 




sp 


where the subscript sp denotes the short period approximation. 

In Ref. [11] the Control Anticipation Parameter (CAP) was defined as: 


CAP 


9o 




(30) 


The CAP is used as an indicator of aircraft maneuverability. Recalling (27) and (25) we get: 
9o 




stick 


UK„ 


stick g 

From equation (24) we can see that a i4 = — and combined with (25) we get: 


(31) 


up q up 

CAP = ^ - — 2E- 


v*. («f)„ 

In Ref. [1] the following limits are defined for level 1 performance: 


(32) 



> 3.5 





0.28 

< 

CAP 

< 

3.6 


V^28f 

< 

Usp 

< 


(33) 

0.5 

< 

Csp 

< 

1.35 



and for (^) 4S < 3.5 
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zH 


h— u 


Figure 2: Basic II ^ Formulation 


( 


< iO 


w sp 


< 


0.5 < Op < 1-35 


(34) 


Based on recent analyses of flight simulations and tests these boundaries have been moved 
closer to one another as follows [12]: 


For 


> 9.0 





\ 0 36 

< 

CAP 

< 

1.2 

< 

lOM 

si < 

a — 

u sp 

< 



[ 0.5 

< 

Cap 

< 

1.1 

and For (• 

Hi') < 

Of ^ 

9.0 



J 

\ 1.77 

< U sp 

< 

3.2 


1 

[ 0.5 

ft, 

VI 

< 

1.1 



(35) 


(36) 


After more evaluations of handling qualities at higher angles of attack [13] larger values of 
the damping coefficient were found to give the best handling qualities, i.e : (, p > 1, and 
a ) sp is recommended to be between the lower bounds of the previous two references. Those 
criteria for and u sp will be used in our definition of level 1 handling qualities. 


4 Hoc Design 

In order to synthesize a controller, the system under consideration is formulated as shown 
in Figure 2, where w is the disturbance vector, z is the error vector, y is the measurement 
vector provided to the controller, and u is the controller output vector. The vectors w and 
z are scaled so that their expected maximum absolute values are less then one. When the 
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controller makes ||T 2 U i||oo < 1, the system is robust with respect to stability and the speci- 
fied performance. With synthesis, only the inputs/outputs related to the performance 
specifications are taken into account. 

The aircraft has two pitch moment control effectors. These are the stabilizers which have 
a deflection limits of -24 to +10.5 Deg and a maximum angular rate of 3 40 Deg/Sec, 
and the Thrust Vector Control (TVC) vanes with deflection limits of -17.5 to +15 Deg and 
deflection rate limits of 100 Deg/Sec. 

The controller output is a moment command that has to be divided between both effectors. 
The deflection command for each of the actuators is computed through a constrained opti- 
mization procedure [14]. Defining the total moment commanded by the controller to be T 
and the total moment produced by the control effectors to be Ms tvc 8 tvc + Ms,8 3 , it is clear 
that 

T = M Stvc S tvc + M s ,6 s (37) 

with M Stvc and Ms, being the aerodynamic moment derivatives with respect to the TVC 
and the stabilizer deflections respectively. We want to allocate these deflections to minimize 
the sum of the squares of the normalized deflections. The cost function for this constrained 
optimization problem is: 




-)’ + ( 77 %-)’ 1 
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tvc\max 


(38) 


Using the Lagrange Multiplier method, we define the extended cost function to be : 


J — J + A(T — Ms tvc 8tvc — M$,8 3 ) 


(39) 


Taking the derivatives with respect to 8 3 , 8 tvc and X and solving for 8 tvc and 8, the resulting 


deflections are: 


T/Ms,,, 

J _| Mh l^lmar 

^ M itve \8tvc\max 


K tvc T 


6s 



K S T 


( 40 ) 

( 41 ) 
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Since the TVC vanes can suffer thermal damage from the jet exhaust if they are deflected for 
long periods of time, there is a wash-out filter with a time constant of 1 Second in front of the 
TVC. The actuator dynamics are approximated by first order systems with time constants 
of 1/20 second for the TVC actuators and 1/30 second for the stabilizer actuators. 

Aircraft Linear Model 

As an example of our design approach, we will design a controller for the flight condition 
of 15000 ft, 20° angle of attack and Mach number of 0.26. The aircraft state vector for the 
design is x = [u, w , q, 6\ , the measurement vector is y = [a, n 2 , <?] and the control 
effector vector is u = [<5, , <5 t „ c ]\ u and w are the body longitudinal and heave velocities 
respectively in ft/sec, q is the body pitch rate in rad/sec, 6 and a are the pitch angle and 
the angle of attack in radians. The normal acceleration in g’s is denoted by n z . 

The linear model of the aircraft at the design point is : 



The open loop system is stable with u p h = 0.136 rad/sec, u ap = 0.620 rad/sec, ( p h = 0.019, 
and ( sp = 0.388. The subscript ph refers to the phugoid mode. 

From (42) we see that Ms, = 5(3,1) = —0.735, Ms tvc = 5(3,2) = —1.683 and Z a = 
A( 2,2) = -0.312. 

The detailed interconnection model for the Hcc synthesis is shown in Figure 3. This figure 
is an extension of Figure 2. The input and output vectors are : 


Z — [ Stvc &tvc *7 Hz A ] 

w = [ stick c n a n n x n q ] 

y = [ stick c q m n Zm a m ]' 

u = T 


(43) 


The first two components of z are the scaled TVC deflection and deflection rate and the 
next two components are the scaled stabilizer deflection and deflection rate. The last three 
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Figure 3: Interconnection Structure for Hoo Synthesis 
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components are the scaled errors in q , n z and a respectively. The w components are the 
scaled stick input and the additive noise in the a , n z and q measurements. The vector y 
is composed of the actual stick deflection and sensor outputs. The vector u is the controller 
output. 

The system is formulated so that all three handling quality criteria can be included in the syn- 
thesis. In this example we only optimize with respect to the angular rate, q. That is we want 
the angular rate to track the pilot stick input. For this case the weighting functions, W Hz 
and W a are set to very small numbers so they will not effect the design objective of sa- 
tisfying the q handling quality criterion. From equations (27-29), the handling qualities 
parameters: u> sp ,K q ,K nz and K a are related to the aircraft linearized model. According to 
the relationship given in (26) (“■”*■)« = 2.5 and from section 3 we choose u sp and ( , p to be 
1 rad/sec and 1 respectively. Our basic requirement is for the aircraft to have the capability 
of performing an acceleration of 1 g due to 1 inch of stick input. This requirement results in 
K nz in equation (28) to be 1 g/inch. Using (25) for w = 0 K q = K nz fj = 0.115 inch . 
From the second equation in (24), for a = 0 we get K a = ^ = 0.368 rad/inch. 

A point to remember is that the handling quality transfer functions are related to the short 
period approximation of the aircraft. For the actual aircraft which includes the phugoid 
mode as well as the short period, the actual steady state values of q and n z for a given stick 
input are zero. 

The inputs to the weighting functions: W q ,W nz and W a are the differences between the 
desired output from the handling quality model, and the actual output from the aircraft 
model. We would like the actual transfer function between the stick and pitch rate to be as 
close as possible to the handling quality transfer function between the stick and desired pitch 
rate, defined in (27). The weighting which defines how well we want to follow the handling 
quality transfer function is W q . We define W q to be: 

w, = 57.3iilfe (44) 

1 ' 3 

This weighting function defines an allowed tracking error of 1 Deg/Sec up to the bandwidth 
of 3 Rad/Sec. Above 3 Rad/Sec, the error is allowed to grow by allowing W q to decrease 
until the frequency of 100 Rad/Sec. We allow larger errors in this frequency range since we 
are not interested in tracking inputs at frequencies greater then 3 Rad/Sec. The weighting 
function does not decrease above 100 Rad/Sec in order to limit high gains at high frequency 
so that high frequency elastic modes effects are not excited. Since we put no constraints on 
n z and a errors, we set W nz and W a to be very small, so the design will not be affected by 
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n z and a handling quality requirements. 


Since n z and a handling quality requirements are derived from the same short period ap- 
proximation as is the q handling quality requirement, we would expect that as long as the 
short period equations approximation is valid, a control law which satisfies the q handling 
quality requirement also does a reasonable job of satisfying the n z and a handling quality 
requirements. This is the case at lower values of angle of attack where there is good separa- 
tion between phugoid and short period modes. As the angle of attack increases, the phugoid 
and short period frequencies approach one another and satisfaction of one handling quality 
requirement does not mean that the others are also satisfied. 

The noise assumed through the design is a white noise with zero mean and with the following 
spectral density components : 0.01 Rad/p^ in a measurement, 0.01 g/yj^~ in n z measu- 
rement and 0.01 0 measurement. The maximum stick deflection for which the 

design was performed is 2.5 inches , the stabilizers and TVC maximum deflections, |6 4 | mai 
and |<S tt)C | max , are 10.5 Deg and 16.0 Deg , and their maximum rates, |<5 a | max and |6<u C | max , are 
40 Deg/Sec and 100 Deg/Sec respectively. Recalling equations (40) and (41) together with 
the aerodynamic coefficients and the above maximum deflections we get that K tvc = —0.500 
and K 3 = -0.215. 

Since the synthesis is a frequency domain design technique, the stick deflection should 
also be characterized in the frequency domain. It is assumed that it takes the pilot at least 
2 Seconds to deflect the stick from zero to its maximum value of +2.5 inches back through 
zero to its minimum deflection of -2.5 inches and then back to zero. This corresponds to a 
bandwidth of 3.14 Rad/Sec. Since the force which the pilot applies is proportional to the stick 
acceleration, the transfer function from the pilot force input to the actual stick deflection is 
a second order system with a bandwidth of 3.14 Rad/Sec. This transfer function is realized 
by the pre — filter block in Figure 3. The transfer function defined for pre — filter is : 


pre — filter 


9.86 

s 2 + 6.28s + 9.86 


(45) 


The above models were used in our initial control law design. We first designed our control 
laws using a short period approximation of the aircraft dynamics. This resulted in a small 
value for 7, but when these control laws were applied to the full aircraft dynamics which 
included the phugoid, very poor performance resulted. Next we decided to include a phugoid 
model in our handling quality requirements. This allowed us to keep 7 small and provided 
excellent performance when the control laws were used with the full aircraft model. In 
addition, if the handling quality model only contains the aircraft short period approximation, 
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Figure 4: -fi-r transfer functions 


the resulting control law design cannot directly modify the phugoid mode. With our approach 
we can stabilize unstable phugoid poles and increase their damping if needed. 

The inclusion of a phugoid mode in the handling quality requirement is done by curve fitting 
the aircraft phugoid approximation at low frequency to the short period handling quality 
specification. Because of the wash out filter in front of the TVC input, only the stabilizers 
are effective at frequencies up to the phugoid frequency. Allowing 10.5 Deg of stabilizer 
deflection for maximum stick deflection of 2.5 Inches we get the low frequency relationship: 


l 

stick 


0.072 


(46) 


The curve fitting and modification to the handling quality specification is shown in Figure 4. 
Note the increased damping of the phugoid mode in the handling quality model compared 
with the actual aircraft transfer function. The maximum singular value over frequency of the 
closed loop transfer function of the system given in Figure 3 is shown in Figure 5. The fact 
that the maximum singular value is less then one means that the system performs according 
to the specification in the different weighting functions that were defined in the Hoc synthesis 
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Figure 5: Closed loop maximum singular value 


formulation. 

Further analyses were done to verify the results. The first was a time simulation of the closed 
loop system applied to the linearized aircraft model. The results are shown in Figure 6. The 
input to the system is a stick doublet shown in Figure 6h. The stick amplitude is 2.5 inches. 
In the figures [variable_hq] stands for the response of the handling quality model to the stick 
input, and [d variable] is the difference between the response of the closed loop linearized 
aircraft model and the handling quality response. Figures 6a and 6b indicate that the 
actuators are almost at their limits, and Figures 6c and 6g show that the actual pitch attitude 
and rate are extremely close to those defined by the handling quality transfer function. Also 
Figures 6e and 6f show that we follow the desired handling quality specifications for a and 
n z even though those were given very small weightings in our design procedure. Figure 
6d indicates that we do not follow the flight path angle from the handling quality model very 
well initially, however the flight path angle is a long period variable whereas our objective 
was to improve short period response. Also as time increases, the flight path angles from 
the handling quality specification and the aircraft axe essentially the same. Thus based on 
linear models our design is successful. 
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Figure 6: Time response of the linear simulation with the H controller 
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Figure 7: Aircraft trajectory on X - Z plane 

Next step in analyzing the Hoc controller is by inserting it into the full nonlinear simulation 
of the HARV. The doublet stick input is identical to that for the linear simulation. The 
aircraft trajectory in x-z plane is shown in Figure 7 and its Mach number is presented in 
Figure 8. The ordinate in Figure 7 shows the deviation from the initial altitude of 15000 ft. 
The Mach number is not continuous because of the way its computation is performed as a 
function of the angle of attack in the simulation. 

It can be seen that the aircraft climbs slightly and then loses about 1000 ft at altitude. It 
then starts to climb again. Figure 8 indicates a decrease in Mach number with the initial 
climb followed by an increase as the aircraft descends. Figures 7 and 8 illustrate the classical 
damped phugoid mode. 

Figure 9 includes the remainder of the aircraft responses. Figure 9a shows the deflection of 
the stabilizers and the TVC vanes. Both reach their maximum values during the maneuver. 
At this simulation no limit was enforced on the TVC vanes so they were slightly above the 
actual limit. Figure 9b shows the stabilizer deflection rate. The TVC deflection rate is 
not an output of the simulation so it is not shown. The deflection rates do not reach their 
maximum values. Thus the maneuver is limited by the control effector deflections rather 
then deflection rates. Figure 9c shows the pitch angle, 6 , from the nonlinear simulation and 
from the handling qualities models, (note that both 6 and a from the handling qualities 
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Figure 8: Mach number as a function of time 

model have been increased by steady state values of 20 Deg which are the equilibrium values 
of these variables for this flight condition). The pitch angle from the nonlinear simulation 
tracks the pitch angle from the handling quality model very closely even though angular 
excursion of approximately 40 Deg from the equilibrium are encountered. The flight path 
angle from the handling quality model is shown in Figure 9d. As in the linear case the 
handling quality flight path angle is not tracked very well, however, tracking of flight path 
angle was not one of our design objectives. Figure 9e shows excellent tracking of angle of 
attack. Normal acceleration is tracked less well. This could be because there is lateral motion 
in the nonlinear simulation which effects the measured normal acceleration. The handling 
quality model is based on the assumption that there is no lateral motion. Figure 9g shows 
that the pitch rate from the nonlinear simulation tracks the handling quality model very well 
and Figure 9h shows that the pitch rate also tracks the stick input profile well. Since the 
control law was designed to provide pitch rate tracking of the stick input, the H 0 0 control 
meets this design criterion very well. 

If robustness issues are ignored we see that the H , <*> design gives us good performance and 
the nonlinear results are very similar to the linear results for relatively large angle of attack 
changes. The H <*> design is very conservative in terms of stability robustness since error 
models are not included explicitly in the synthesis procedure. The \i synthesis procedure 
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includes error models and is the subject of the next section. 


5 fi synthesis controller 

The n synthesis controller was designed to compensate for the following four sources of 
uncertainty in the math model of the aircraft: 

1. Parametric uncertainties in the linear aircraft A matrix. 

2. Actuator effectiveness uncertainty. 

3. Measurement uncertainty. 

4. Input uncertainty. 


The first uncertainty will be presented as an LFT while the others will be treated as mul- 
tiplicative uncertainties. All the uncertainties will be translated into an extended LFT in 
order to perform the fj, synthesis. The synthesis software available at this time assumes all 
the perturbations and uncertainties are complex. A real fi analysis will be presented along 
with the complex fi analysis and we shall see that both axe very close to each other. This 
indicates that very little performance is lost by modeling uncertainties as complex functions. 

Since our controller will not be scheduled on any of the measurements, we shall carry out 
our design for a small range of Mach number and angle of attack variations. 

The three points in the flight envelope for which the controller will be synthesized are: 1) H 
= 15000 ft, Mach = 0.291, a = 15° 2) H = 15000 ft, Mach = 0.258, a = 20° and 3) H = 
15000 ft, Mach = 0.246, a = 25° 

The linear model at 20° is given by equation (42) while the other two linear models are given 
below. 


At H = 15000, Mach = 0.291 and a = 15° the linear model is: 


A15 

B 15 

Cl5 

D 15 


0.004 

-0.132 

0.0004 

0.0 

- 0.001 

0.004 

0.0 


0.025 

-0.315 

- 0.0021 

0.0 

0.003 

0.010 

0.0 


-78.939 

294.589 

-0.199 

1.0 

0.0 

0.091 

1.0 


-31.071 

-8.327 

0.0 

0.0 

0.0 

0.0 

0.0 


-1.271 

-12.033 

-0.862 

0.0 

0.0 

0.376 

0.0 


2.771 

-15.116 

- 1.866 

0.0 

0.0 

0.470 

0.0 


( 47 ) 
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Figure 9: Time response of the nonlinear simulation with the Hoo controller 
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Figure 10: A22 variations as an LFT 


At H = 15000, Mach = 0.246 and a = 25° the linear model is: 


^25 

5° 

Cn 

C 25 

D 25 


0.024 

0.002 

-109.202 

-0.116 

-0.286 

234.186 

0.0002 

- 0.0001 

-0.167 

0.0 

0.0 

1.0 

- 0.002 

0.004 

0.0 

0.004 

0.009 

0.070 

0.0 

0.0 

1.0 


-29.153 

-1.520 

0.432 

-13.596 

-7.529 

-11.188 

0.0 

-0.604 

-1.367 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.239 

0.348 

0.0 

0.0 

0.0 


(48) 


The elements, A22 A 23 A32 and / 1 33 , in each of the above linear models define the short 

period mode. The nominal system for the design is for the 20° flight condition. 

As an example of modeling the parametric variations in A let us examine how A22 variations 
can be brought into LFT form. We shall define A22 nom = —0.312 , A22 min = —0.315 and 
A22 max = —0.286. The variations in A22 are shown in Figure 10 in LFT form. The relation 
of A22 to <5 is: 


A22 — T 0 i 


-0.312 + 


0.0054(5 
1 - 0.793(5 


(49) 


One can check that for <5 = 1 T ot = —0.286, for 6 = 0 T 0 » = —0.312 and for <5 = — 1 T 0 , = 
—0.315. Variations in A 23 , A32 and A33 are modeled in a similar manner. Variations in the 
remaining elements of A are not considered in the design. 
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Figure 11: 6 S ,6 tvc and q m variations. 


Putting the variations of the four elements in A which are related to the short period mode 
into an LFT form gives us the following structure: 



Where Z\ Wi are related to A22 variations, z 2 w 2 to A23 variations, z 3 w 3 to A 32 variations 
and 24 W4 to A 33 variations. Every z t is connected to its corresponding w, through |<5, [ < 1 
where i = 1,2, 3, 4. 

Assuming a 25% multiplicative uncertainty in the effectiveness of the control effectors and a 
5% multiplicative uncertainty in the measurement scaling, the blocks that should be inserted 
into Figure 3 are shown in Figure 11. The 6 a and 8 tvc blocks should be inserted at the aircraft 
dynamics input, and q m at the controller input. 

Input uncertainty was used for modeling the uncertainty of the system dynamics at both 
low and high frequency. The input uncertainty model is based on the assumptions that 
1) the uncertainty at low frequencies is 20% 2) the first elastic mode is at 30 Rad/Sec and 
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Figure 12: Input uncertainty structure. 


3) the elastic modes add 600% of uncertainty to the system dynamics. The resulting input 
uncertainty model was selected as: 


input uncertainty 


6(3 + 1 ) 
s + 30 


(51) 


With this uncertainty model, the structure shown in Figure 12 has to be inserted after the 
controller output in Figure 3. Taking all the above uncertainties and perturbations into 
consideration, we define a new system interconnection in the format of Figure 2. In this case 
y and u are the same as those defined in equation (43) for the H^ synthesis. The z and w 
vectors for the p synthesis are: 


z = [zi z 2 Z 3 z 4 z 5 z 6 z 7 z 8 S tvc 6 tvc S 3 S s q h z a ] 

w — [ujj w 2 W 3 W 4 u> 5 w 6 W 7 Wg stick c n a n„ r h q ] 


(52) 

(53) 


Note the difference in the dimensions of w and z compared to those defined in equation (43) 
for the Hoo synthesis. The additional inputs/outputs are necessary for addressing stability 
robustness. 

Using the above definitions for all the inputs and outputs of the basic interconnection struc- 
ture shown in Figure 2, we went through thq D-K iterations. Starting with the first 7 of the 
Hoo design to be 25 and the first p value of the p analysis to be 10 we went through 7 D-K 
iterations to get the final and lowest p we could achieve. 

The final p achieved is shown as a function of frequency in Figure 13. The full line is the 
upper bound for the complex p analysis. Its maximum is 2.77. The dashed line is the upper 
bound for the real p analysis, and its maximum is 2.72. 

The p controller was applied to the three linear models of the aircraft at a = 15° , a = 20° 
and a = 25°. Only one unstable pole resulted. This was for a = 25° and it had a very slow 
time constant of 1250 Sec. Using the Hoo controller with the same linear models we get an 
unstable pole with a time constant of about 100 Sec at a = 25°. 
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Figure 13: complex and real /z analyses results. 
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To test for the ^ controller performance, 6 linear simulations were done. 3 with the // 
controller and 3 with the controller. Each controller was simulated with the aircraft 
linear models at 15°, 20° and 25°. The results of comparing the pitch rate, q, in all the 
cases are presented in Figure 14. The only point the Hoo controller has a better response is 
at the nominal design condition at a = 20° . At the other two points the fx controller is 
better then the H <*> controller in the sense that its settling time is smaller, its response is 
less oscillatory and it can better track the handling quality response. 

The \x controller that was synthesized has 111 states. Using a system balanced model reduc- 
tion we designed a 25 state controller with the upper bound on // being 2.82. This reduced 
order controller was inserted into the full nonlinear HARV simulation. The nonlinear results 
along with the handling qualities responses are given in Figure 15. 

Even though no major differences can be noticed between both nonlinear simulations — the 
one with the H <*> controller, Figure 9, and the other with the fx controller, Figure 15 - the ji 
controller is superior to the Hoo controller from the linear analyses. As mentioned before, the 
closed loop system with the fi controller has a much slower unstable mode then the closed 
loop system with the Hoo controller. The system with the fx controller performs better at all 
points other than the nominal design point. 

In the above we showed how we succeeded in synthesizing a // controller that addresses the 
stability robustness and at the same time provides performance close to the specifications 
with the three linear models as well as with the nonlinear simulation. When adding a 
scheduling algorithm on a or on the Mach number, the maximum value of fx should be lower 
then 2.77. Hopefully this will result in fx close to one. 

Adding the variations in the other elements of A matrices to capture the phugoid modes 
variations will make the D-K iteration longer. This is not necessary as long as the phugoid 
modes are stable, as they are in our case. 


6 Controllers Comparison 

Before implementing a dynamic inversion control law on the HARV, we are synthesizing a 
fi controller that will follow the a handling quality transfer function as defined in equation 
(29). The angle of attack handling quality transfer function will be followed at the high 
angles of attack regime. 
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Figure 15: Time response of the nonlinear simulation with fi controller 
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Before synthesizing the fi controller we synthesized an H <*> controller that should meet only 
the performance specifications at a = 35 Deg. In this section we shall compare our controller 
(UM) performance with the baseline controller performance in a maneuver that includes a 
full throttle and a longitudinal stick input as described in Figure 16i. A much higher stick 
deflection is needed for trimming the aircraft with the baseline controller relative to UM 
controller. Since the aircraft goes through large excursion in angle of attack and Mach num- 
ber the nonlinearities in the equations of motion play an important role in these simulation. 
Also as shown in Figure 161 there is a substantial roll angle during the maneuver, so the 
lateral dynamics also are important. 

Initial the aircraft is in trim at a = 35 Deg , H = 15000 ft and at Mach = 0.205. As can be 
seen from Figure 16 both controllers are very close to each other. 

This example is one of many comparisons we did between our control laws and the baseline 
control law. Our controllers performed as well as or better then the baseline controller in all 
cases. 
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Figure 16: Baseline Controller(full) and UM Controller(dashed) 



















7 Conclusions and Plans for Further Research 


The methodology outlined in this report has wide application to design of H 0 c based control 
laws for aircraft. Both the H 0 Q and // synthesis control laws provide excellent dynamic re- 
sponse characteristics. Once the basic techniques have been developed control laws can be 
designed quickly and these control laws provide dynamic response characteristics which are 
as good or better than those of the baseline control law. The p synthesis control laws show 
better performance at off nominal flight conditions than do the H 0 c control laws when eva- 
luated using the linear simulation, but both control laws yield essentially the same response 
in the nonlinear simulation. Addition work is needed in evaluating the robustness properties 
of these two control laws. We also need to study the transitioning of the control laws from 
a pitch rate command to an angle of attack command as angle of attack increases. We are 
beginning to work on the design of the nonlinear dynamic inversion control law. 
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